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1. Introduction 

The axial charge of the nucleon, or more precisely the ratio Ga (q 2 = 0)/G v (q = 0) has been 
determined to a high precision from neutron j8 decay, with Ga (0) / G v (0) = 1 .2695 (29) . In general, 
the axial form factor Ga.bb' f° r an octet baryon is given by 



where Gp is the induced pseudoscalar form factor. The axial charge is defined as the value of the 
axial form factor at zero momentum transfer GA,BB'(q 2 = 0). m the following, we will omit the 
indices B and B' when referring to the nucleon. For the nucleon in the chiral limit, the Goldberger- 
Treiman relation Ga = PngnNN /M^ connects the axial charge to the pion decay constant F n , the 
pion-nucleon coupling constant g^N an d the nucleon mass M#. Away from the chiral limit, this 
relation is still approximately fulfilled. Assuming the conservation of the vector current (which 
is the case for mass-degenerate light quarks m u = m<i), the nucleon axial charge is also related to 
the polarized quark distributions in the proton: Ga = Aw — Ad In an isovector combination, 
disconnected contributions cancel, making high-precision lattice computations feasible. 

The Chiral Perturbation Theory {%PT) expressions relevant to the nucleon axial charge have 
been calculated in [||], where finite volume effects are taken into consideration. While a recent 
simulation with domain wall fermions |^] finds considerable finite volume effects and scaling in 
M K L, volume effects calculated in %PT lead to differing conclusions. Trying to attribute this dif- 
ference to excited state contaminations arising from finite separation in Euclidean time, Tiburzi [Q] 
estimates the effects of such contaminations and obtains that they would lead to an over-estimation 
of Ga rather than an under-estimation. He also suggests to study Ga using the variational method. 
Lattice results for the nucleon axial charge have furthermore been presented in For a recent 
review, please refer to the review by Renner [0j. 

So far, only one group has reported results for the axial couplings of sigma and cascade hy- 
perons [Q]. The corresponding #PT calculations can be found in [§]. In [§] input from experiment 
and lattice QCD is used to determine the unknown parameters in the %PT expansion and predict 
the mass dependence and values of the axial charges in the chiral limit. 

In the next section, we explain the setup for calculations of baryon axial charges using Chirally 
Improved (CI) lattice fermions and the variational method. We will then move on and present 
results from our calculations of the axial charges of the nucleon and of £ and E hyperons. 

2. Details of our calculational setup 

For our simulations we use CI fermions [|10|], which are approximate Ginsparg -Wilson fermions 
based on an expansion of Dirac operator terms on a hypercube. CI fermions have been tested ex- 
tensively in quenched calculations [Jll]] and results for the ground state spectrum of mesons and 
baryons from dynamical CI simulations have been presented recently [|L^]. 

Assuming mass-degenerate up- and down quarks it is sufficient to consider the following cur- 
rent insertion to extract the axial charge 




(1.1) 




uy^y 5 u 



(2.1) 
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t = t' t" t = t\ < t' < t 2 t" t = t' t" 

Figure 1: Different possibilities for the calculation of sequential quark propagators (left and middle). The 
right-hand side shows an illustration of the full baryon three-point function. 



in which u denotes an up quark and d denotes a down quark. In the following, we will show how 



three-point functions with insertions like those in Equation [2JJ can be evaluated on the lattice. 

For our calculations we use so-called sequential propagators. In [|IJ], two methods for the 
calculation of sequential quark propagators are presented. Figure [j] illustrates these two approaches. 
On the left-hand side, sequential sources are built from specific diquark propagators, while in 
the middle, the propagators are calculated for each possible insertion separately. An illustration 
of the full baryon three-point function is provided on the right-hand side. Which approach is 
computationally cheaper depends on the physics objective. In our case, we want to use two different 
insertions and two widths of smearing for three different interpolator types. As we are using a 
rather coarse lattice, a small number of insertion timeslices should be enough. Therefore, even 
when just considering the nucleon case, the second approach using sequential sources from single 
quark propagators is slightly cheaper. Moreover, these propagators can subsequently be used for 
other hadrons and for the calculation of transition form factors. 

To relate lattice operators, which receive a finite renormalization, to their continuum coun- 
terparts, we need to estimate the renormalization factors Zr of the bilinear currents in question. 
In general, we have to multiply the lattice result Gj? by the appropriate renormalization factor to 
obtain values that can be compared with results extracted from experiments 

G p T hys =Z r G^, (2.2) 

which are typically given in the modified minimal subtraction (MS) renormalization scheme. 

For dynamical CI fermions, these renormalization constants have been estimated using local 



bilinear quark field operators in [ 13p. It would however be useful to have an independent estimation 



of these constants from a different method. In the case of the vector current, one can estimate the 



constant Zy by calculating the vector charge Gy denned in analogy with ([LI]) via 

(B'\V^q)\B) = u B ,{p') ^ Gv {q 2 )+q v a vtl ^pj u B {p)^ x , (2.3) 

as Gy(q 2 = 0). This quantity has to be 1 in the continuum, as it is related to the electric charge of 
the proton in the limit of equal quark masses |JlJ|. 

For lattice fermions with exact chiral symmetry, the axial vector renormalization constant Za 
and the vector renormalization constant Zy have to be equal. For lattice fermions which only fulfill 
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the Ginsparg-Wilson relation approximately, there should be small deviations from this. To obtain 
an independent estimate of Zy, we use a ratio of two-point over three-point functions [14] 



R (k) 



]J m 



(*) 



(2.4) 



where C(t) is the matrix of two-point correlation functions and Ty(t,t') is the matrix of three-point 
correlators with a vector insertion. The eigenvectors y/ are the ones obtained from a variational 
analysis of C(t). We then compare with the preliminary estimates from [13]. The ensemble names 
are according to [|l2|], where details of the run parameters are provided. For runs A and B the two 
methods agree within 2 — 3%. While a determination using local quark bilinears yields 0.818(2) for 
run A and 0.826(1) for run B, we find values for Zy of 0.803(2) and 0.792(2) respectively. For run 
C there is a rather large discrepancy and the method of [13] leads to a value of 0.829(1) while we 
obtain 0.77(1). Notice also that two different methods for the determination of the renormalization 
constants are presented in [ |l3| ] which only agree after a chiral extrapolation of the results is per- 
formed. At the same time, the ratio Z^/Zy determined from the values in [13] is almost identical 
for both methods used and also stable under chiral extrapolation of the results. 

In our determination of the axial charge from run C, we encounter what we suspect to be large 
finite volume effects. Notice that the value of Zy obtained from the nucleon three-point functions 
might be plagued by the same effects. As we cannot calculate Za from baryon three-point functions, 
we therefore always use the ratio Z^/Zy from [13]. In the next section, we discuss in detail which 
ratios we measure on the lattice to obtain the renormalized axial charge Ga- 



3. Nucleon axial charge from dynamical CI fermions 

The usual approach [JTJ, |3[] is to extract the nucleon axial charge from ratios of Ga over Gy 

Z A T^(t,t') 

GA = Z-yTyW)' 

using single correlation functions built from either smeared quarks or gauge fixed box or wall 
sources. This approach has the advantage that some of the systematic errors entering the lattice 
determination will cancel. 

We instead use the variational method, which is commonly used to extract ground and excited 
state masses. It is based on a correlation matrix Qj(t) = (0,-(f)0j(O)) where are operators 
with the quantum numbers of the state of interest. The eigenvalues A, of the generalized eigenvalue 
problem C(t)vj = A,-C(?o)v,- may be shown to behave as A,-(f ) <x eT tE ' (l + G (e - '^')) , where is 
the energy of the i-th state. The approach may be generalized to three point functions. Following 



Ql4|], we obtain an expression for Ga'. 



^Zt JUZlVrfnMtlVrf (32) 

Zv LlLmWl k) Mt,t')lmY^ ' 

Figure ^| shows a typical plateau for the axial charge of the nucleon from run C extracted from such 
a ratio. The horizontal lines denote the results from a linear fit in the displayed range. Notice that 
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Figure 2: Example plot to illustrate typical plateaus observed with our variational basis for a source-sink 
separation of rs 1 .2fm. The data is from run C. 
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Figure 3: We compare our results for Ga to a recent determination from domain wall fermions. Our data is 
labeled as "2 flavor CI". Results from 2+1 flavor domain wall fermions are taken from Yamazaki et al. [g]. 
On the left-hand side, we plot the results over M\. Towards lower quark masses finite volume effects are 
clearly visible. On the right-hand side we display our data in units of M K L. 



we observe a plateau in the full range of points we calculated. For all three ensembles, we choose 
timeslice t = 9 for the position of the sink. This corresponds to a source-sink separation of roughly 
1.2fm. For run B, we currently only have data for insertion timeslices t' from 5 to 9. Instead of 
assuming that the central value at 5 is the physical one, we perform a linear fit in the range 5 to 7. 

We compare our data to recent results from domain wall fermions ^ in Fig. ^[ The left-hand 
side plot shows the results for Ga plotted over the square of the pion mass M\. While results at 
large pion masses lead to values close to the experimental value, the result from run C deviates 
substantially from this behavior. The same is true for the domain wall data and this behavior seems 
to be a universal feature associated with finite volume effects [Q, On the right-hand side of the 
figure we therefore plot the results for Ga over M K L, where L corresponds to the spatial extent of 
the lattice. This plot can be directly compared to Fig. 3 of 

Before we move on to calculations for hyperons, let us briefly comment on the sink-dependence 
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Figure 4: Results for the axial charges of the £ (left-hand side) and S (right-hand side) hyperons compared 
to the mixed action results by Lin and Orginos M\. 



of our results. While results from run A and B are rather insensitive to the sink location in the region 
explored (timeslices 9-13), a systematic shift upwards can be observed for run C when reducing the 
separation between the source and the sink from 8 to 6 timeslices which corresponds to distances 
of 1.2fm and 0.9 fm. We want to point out that this does not affect the quality of the plateau which 
still stretches over the entire region of insertion times. Taking a look at the nucleon two point func- 
tions, contributions from excited states to the ground state of the variational analysis are visible up 
to timeslice 4. This is an indication that excited states may indeed be responsible for measuring a 
larger value of Ga if excited state contributions are not sufficiently suppressed. With just 50 config- 
urations, the statistical errors from our preliminary dataset are by far too large to make a stronger 
and more quantitative statement. 

4. Hyperon axial charges 

In this section we present results for a calculation of hyperon axial charges. For the £ and E 
hyperons we adopt the following definitions: 

(Z + \Al\L + )-(L-\Al\L-)=Gzzu v y^ 5 u v , (4.1) 
(E \Al l \E°)-(E-\Al\E-) = G EZ u v r,r 5 u v . (4.2) 

Again, no disconnected contributions appear in the isovector quantities and the calculation proceeds 
similar to the nucleon case. In particular, no additional sequential propagators are needed. 

Figure |] shows our results for the axial charges of the £ and S hyperons. We compare our data 
to [f7|] and we can see a quantitative agreement in the full range of masses. Unlike for the nucleon, 
no significant decrease is observed towards the chiral limit in the plot for the £ (l.h.s.). The plot 
for the S (r.h.s.) also shows a nice agreement with the results from In this case our value 
corresponding to the smallest pion mass shows a slight decrease towards the chiral limit, but the 
error bars are large and this may as well be an effect of our limited statistics. Our purely statistical 
errors on the preliminary dataset of 50 configurations are still large but can be substantially reduced 
by using our full statistics. 



6 



Baryon axial charges from CI fermions 



Daniel Mohler 



5. Summary and outlook 

We have presented preliminary results from a calculation of baryon axial charges using a 
full variational basis to efficiently suppress contaminations from excited states. We used a basis 
of baryon interpolators with different Dirac structures and two different smearing widths for the 
quarks. The results are in good agreement with the literature and we obtain clear plateaus for ratios 
calculated with the method of [|14|]. Provided the signal for the states in question is strong enough, 
this method can also be applied to several other quantities of interest, among them the axial charge 
of the delta baryon and the N-A or E-A transition. 

In general, the method we use can also be applied to three-point functions involving excited 
states, provided that the signal is good enough to ensure the necessary separation between the 
source/sink and the current insertion. 
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